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We investigate multiple Fano, comblike scattering resonances in single-layered, concentric core- 
shell nanoparticles composed of magnetic materials. Using the Lorenz-Mie theory, we derive, in 
the long-wavelength limit, an analytical condition for the occurrence of comblike resonances in the 
single scattering by coated spheres. This condition establishes that comblike scattering response 
uniquely depends on material parameters and thickness of the shell, provided that it is magnetic 
and thin compared to the scatterer radius. We also demonstrate that comblike scattering response 
shows up beyond the long-wavelength limit and it is robust against absorption. Since multiple Fano 
resonances are shown to depend explicitly on the magnetic permeability of the shell, we argue that 
both the position and profile of the comblike, morphology-dependent resonances could be externally 
tuned by exploiting the properties of engineered magnetic materials. 
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I. INTRODUCTION 

Optical systems exhibiting comblike resonances have 
a wide range of scientific and technological applications, 
especially in spectroscopy and frequency metrology fj]- 
[§]. Comblike profiles, which consist of ultrasharp res¬ 
onance peaks, can be achieved e.g. in the electromag¬ 
netic (EM) scattering by plasmonic nanoparticles or via 
morphology-dependent resonances in nonmetallic scat- 
terers. Indeed, plasmonic nanoparticles support local¬ 
ized surface plasmon resonances that are strongly de¬ 
pendent on their geometrical and material parameters, 
allowing for engineering the scattering response [J. In 
addition, scatterers with a high degree of symmetry, 
such as spheres, spheroids and cylinders, may support 
morphology-dependent resonances, which are related to 
constructive interferences confined inside the particle by 
almost total internal reflections Q. 

Recently, a great deal of attention has been devoted 
to Fano resonances in plasmonic nanostrutures due to 
its sensitivity to both geometry and local environment 
changes The interference between a broad bright 

resonance and a narrow dark resonance modes, supported 
by plasmonic nanostrutures, gives rise to the Fano reso¬ 
nance, which has a characteristic, asymmetric lineshape 
with a narrow bandwidth )7|]. By properly choosing the 
materials and/or the system design, it is possible to gen¬ 
erate e.g. multiple Fano, comblike scattering resonances 
in plasmonic nanoparticles @-0. The design of these 
Fano-comb nanoparticles may enable several applications 
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at different spectral ranges, such as improving the reso¬ 
lution in comb-spectroscopy techniques Q, and optical 
tagging 0. 

One possible way to achieve comblikc resonances in EM 
scattering is to consider multilayered core-shell nanopar¬ 
ticles 0 ] . The main idea is to combine localized plasmon 
resonances and the scatter ing cancellation technique for a 
single core-shell scatterer 0-0 , whose core is dielectric 
and the coating is composed of multilayered semiconduc¬ 
tor materials with a gradient of doping level. An alter¬ 
native approach to obtain Fano-comb resonances is to 
consider nonconcentric, single-layered cylindrical nanos¬ 
tructures [0, whose fabrication is considerably simpler 
than multilayered nanoparticles. 

In the present paper, we propose an alternative ap¬ 
proach, based on magnetic materials, for generating 
comblikc scattering response in concentric, single-layered 
core-shell nanoparticles, which can exhibit Fano reso¬ 
nances 0-0], Instead of inducing Fano-comb reso¬ 
nances by the geometrical symmetry breaking in single¬ 
layered core-shell nanostructures 0, we suggest the 
use of high permeability covers [0 in center-symmetric 
scatterers. Indeed, in the long-wavelength limit, the 
presence of a ferromagnetic layer breaks the isotropy of 
the Rayleigh scattering, and allow us to obtain multi¬ 
ple Fano, comblike scattering resonances composed of 
both electric and magnetic dipole resonances in the ex¬ 
tinction cross-sections [0. These Fano-like resonances, 
known as “unconventional” Fano resonances (0, are dif¬ 
ferent from the conventional Fano ones [S], since they are 
polarization-independent and arise from the interference 
between EM modes with the same multipole moment. 

We explore these properties to derive analytical condi¬ 
tions for the occurrence of comblike resonances in coated 
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spheres composed of a dielectric core and a high magnetic 
permeability shell. As the condition for the occurrence of 
multiple Fano resonances depends explicitly on the value 
of the magnetic permeability of the shell, we argue that 
both the position and profile of the comblike resonances 
could be externally tuned. This tuning could be achieved 
by exploiting e.g. the properties of engineered magnetic 
materials, where the magnetic permeability can be mod¬ 
ified by applying an external magnetic field and/or by 
varying the temperature. 

This paper is organized as follows. In Sec. [Til 
we briefly review the Lorenz-Mie theory for center- 
symmetric coated spheres [2(| [27} ■ By applying the long- 
wavelength approximation, we obtain in Sec. IIIII an ana¬ 
lytical condition to determine the position of the multiple 
Fano, comblike scattering resonances as a function of the 
shell parameters. In Sec. IIVI we go beyond the long- 
wavelength limit and consider a system composed of a 
dielectric nanosphere coated with a magnetic thin shell 
whose permittivity is provided by a realistic, lossy Drude 
model. Finally, in Sec. E we summarize our main results 
and conclude. 


II. THEORETICAL BACKGROUND: 

LORENZ-MIE FORMALISM 

Here we present the main results of the analytic so¬ 
lution of EM scattering by coated spheres, kwown as 
the Lorenz-Mie theory (26|, that are extensively used 
throughout this paper. Let us assume for the incident 
EM wave the time harmonic dependence e~ lut , with ui 
being the angular frequency. The scatterer is a center- 
symmetric, core-shell sphere, with inner radius a and 
outer radius b , composed of linear, spatially homoge¬ 
neous, and isotropic magneto-dielectric materials. The 
electric permittivities (e) and magnetic permeabilities 
(fi) of core (0 < r < a), shell (a < r < b) and embedding 
medium (b < r ) are (ei,/Lti), (£ 2 ,p 2 ) and (£o,/xo)j respec¬ 
tively. The scatterer geometry is depicted in Fig. [TJ 



sidered to be the vacuum and the relative parameters are 
e qr = £ q /eo and y ql = y q /with q = 1 for the core and 
q = 2 for the shell. The scattering efficiency (which is 
the respective cross-section in units of nb 2 ) is 

r, 00 

Q SC a=-^(2n+l)(|a rl | 2 + |& rl | 2 ) , (1) 

y n=i 


where y = kb is the size parameter of the outer sphere 
(k = 2tt/X being the incident wave number) [2(|. The 
extinction and absorption efficiencies are, respectively, 
Qext = (2/y 2 ) E^i( 2n + l)Re(a n + b n ) and Q abs = 
Qext — Qsca- Here, the Lorenz-Mie scattering coefficients 
a n (electric) and b n (magnetic) are [2(| [ 27 } : 


(D n /m 2 +n/y)tpn{y) - ^n-i{y) 

a n = ~ V i (2) 

yD n /m 2 + n/yj £ n (y) - £„-i (y) 

( fh 2 G n + n/y) 4> n (y) - V'n-i (y) 

bn = ^—3 - 4 - , (3) 

yn 2 G n + n/yj £ n (y) - £„_i (y) 
where one defines the auxiliary functions [26} 


Pn{m 2 y) - A n x' n (m 2 y)/ip n (m 2 y) 

1 ^ A n Xn(rn 2 y)/4> n (m 2 y) 

Pn{m 2 y) - B n x' n (m 2 y)/ip n (m 2 y) 

1 - B n Xn(m 2 y)/ip n (m 2 y) 

ip n (rn 2 x) [m 2 P n [miX) - fh 1 P n {m 2 x)} 
m 2 P n (mix)xn{m 2 x) - fhix' n (m 2 x) 
'ip n (m 2 x) [fh 2 P n (m 2 x) - fhiD n (mix)] 
fh 2 x' n {m 2 x) - fh 1 D n (mix)xn{m 2 x) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


with x = ka being the size parameter of the inner sphere 
and P n {p) = d[ln^> n (/9)]/dp. The functions ipn{p) = 
pjn(p), Xn(p) = ~py n {p) and £ n (p) = ijj n (p) - iXn(p) are 
the Riccati-Bessel, Riccati-Neumann and Riccati-Hankel 
functions, respectively, with j n and y n being the spherical 
Bessel and Neumann functions [2(|. The relative refrac¬ 
tive index is m q = y/£ qr p qr and fh q = ^£ ql /p qT , with 
q = { 1 , 2} [ 23 }. Notice that fh q is the inverse of the usual 
definition of relative impedance, Z ql = y, / u„ r /£„/, so that 
fh q = m q for nonmagnetic materials (p ql = 1) [23}. 

Equations © © are the exact expressions for core¬ 
shell spherical scatterers with arbitrary geometrical and 
material parameters. In the following, we consider ap¬ 
proximations for both core and shell parameters regard¬ 
ing the incident wavelength. 


Figure 1. (Color online) A core-shell scatterer with spherical 
symmetry irradiated by a plane wave. The core has optical 
properties (ei, fii) and radius a, whereas the shell has (£ 2 , P 2 ) 
and radius b. The surrounding medium is the vacuum (£ 0 , yo). 

For the sake of simplicity, the medium (£q,Mo) is con- 


III. ANALYTICAL RESULTS 

We proceed our analysis by taking into account the 
long-wavelength regime (y = kb <C 1). Within this ap¬ 
proximation, the arguments m\x (core) and m 2 y (shell) 
of spherical Bessel and Neumann functions should not 
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necessarily be small. Following Tribelsky et al. [25| . 
one can write a n = /[J 7 ^ 1 + tGn^] and b n = 

J-n 3 ' 1 /[J 7 n B ' > +iSn B ' > }- For coated spheres, we calculate 
the new auxiliary functions in the small-particle limit; 
they read 


rl A) 


S ( n A) 


(2n + l)!! L 
(2n- 1)!! 
y n 
y n 

(2n + l)!! L 

yU+2 


(n + 1 )m 2 - yD r 
rh 2 n 


D 


y 


2(2n + 3)!! L 


m 2 yG n - (n + 1) 


m 2 yG n - (n + 3) 


Si B) « - 


(2n — 1)!! 


m 2 G n 


( 8 ) 

(9) 


( 10 ) 

( 11 ) 


We find that there are electric and magnetic scatter¬ 
ing resonances (i.e., constructive interferences) in \a n \ 2 
and |& n | 2 when Qi A ’ B ' 1 = 0 (|a n | 2 = \b n \ 2 = 1), respec¬ 
tively. Similarly, destructive interferences occur when¬ 
ever J 7 L A,B ' > = 0 (|a„| 2 = \b n \ 2 = 0). These destructive 
interferences lead to very small scattering cross-sections, 
and characterize an EM “cloaking” mode. Note that the 
relevant variables to analyze electric and magnetic res¬ 
onances are indeed Gn^ / and Gn^/J 7 ^, respec¬ 
tively §■ For nonmagnetic particles (y i r = y 2l = 1) 
the leading term in Eq. m vanishes, so that it requires 
higher accuracy relative to other approximatio ns [258 . In 
particular, the expressions presented in Ref. [251 ] for a 
homogeneous sphere are retrieved as one performs the 
following mapping: a —> b [D n = D n (m\x)m 2 / m\, 
G n = D n (rnix)fhi/m 2 ] or a — > 0 or m\ = m 2 [D n = 
G n = D n {m 2 y)). 


A. Shell in the ferromagnetic limit (\m 2 \y 1, 

V < 1) 

In addition to the small-particle limit (y -C 1), if 
|m 2 1y 1 (i.e., large permittivity and/or large perme¬ 
ability shell materials), one has ipn('fn 2 y) Ki sin(m 2 7/ — 
nn/2) and Xri(rn 2 y) ~ cos(m 2 y — mt/2) [26|. To simplify 
our analysis, let us define the “phase shifts” 

Ai A) = tan _1 (—A„) , (12) 

Ai B >=tan (13) 

where A n and B n are defined in Eqs. © and ( 0 ), respec¬ 
tively. The ferromagnetic limit (\m 2 \y 1) is obtained 

from Eqs. © ED by replacing D n with cot[777 2 2/+A(( 4 ) — 

-—- / D\ 

7i7r/2] and G n with cot [m 2 y+ A/ — mr/2], and redefining 
the auxiliary functions. In particular, the dipole approxi¬ 
mation (n = 1) for the ferromagnetic small-particle limit 


leads to: 


V (A) ~ V / 

1 ~ 3 l 


Si 


(A) 


•'T 


(B) 


G\ 


( B) 


2 m 2 


1 fm 2 
y\y 

-| |m 2 y + 2 cot 

1 1 

— < m 2 -cot 

V l V 


ytan TO 22 / + A^Jj , 
tan m 2 y + A 1 / 4 ' | , 

m 2 y + A[ B) | , 

™ i A ( b Y 

vn 2 y + A( 


(14) 

(15) 

(16) 
(17) 


In Eqs. (ED ED, the dependence on the core param¬ 
eters (ei,/ri) and x = ka is encoded in A^ A,B \ Re¬ 
garding the shell parameters (£ 2 ,^ 2 ), we consider two 
reciprocal configurations that fulfill the limiting case of a 
shell with high refractive index (|m 2 | = |^/£ 2 rM 2 r| ^ 1), 
namely: (i) large permeability shells with small permit¬ 
tivity (|£ 2 r| ■C l;|/X 2 r| ^ 1) and (ii) large permittiv¬ 
ity shells with small permeability (|£ 2 r| 1; |/^ 2 r| "C 1). 

By symmetry, one case can readily be obtained from the 
other by replacing e with y for our purposes, we focus 
on case (i), i.e., high permeability materials with small 
permittivities such that \y 2r \ O' |£ 2r |- 

We recall that vanishing shell parameters facilitate the 
scattering cancellation (EM cloaking). This mechanism 
is based on the cancellation of the EM fields due to a lo¬ 
cal negative polarizability coefficient. The electric dipole 
(n = 1) is associated with the local polarization vector 
P(r, u>) = £o[£ r (r, u) — l]E(r, cj), where one considers the 
local fields and parameters inside the scatterer. Thereby, 
in case (i), one may locally induce a polarization vector 
out of phase to the electric field, allowing for the par¬ 
tial cancellation of the scattering efficiency (Q sca ~ 0). 
In addition, large values of the magnetic permeability 
in small-particle limit are responsible for breaking the 
isotropy of the Rayleigh scattering, favoring backscat- 
tering [23j] and EM resonances inside the particle [23 ] . 
Consequently, one could in principle engineer suitable 
properties to tune both cloaking and strong scattering 
responses in magnetic core-shell spheres. 

In case (i), spherical shell materials with \y 2l \ 1 

lead to high refractive (|m 2 | ^ 1) and low impedance 
(1 777 . 2 1 1) indices, since low permittivity is consid¬ 

ered. Therefore, from Eqs. ED and ED, the electric 
and magnetic cloaking in the ferromagnetic limit occur 
for t&n[m 2 y + A^f 4 '] ss 0 and cot[m 2 y + A^ S) ] ss 0, re¬ 
spectively. Conversely, for case (ii), |£ 2 r | ^ 1 leads to 
|m 2 | 1 and 1 777 . 2 ! ^ 1, which readily provide the electric 

and magnetic cloaking conditions: cot[m 2 y + A^‘ 4 '] « 0 
and tan[m 2 y + A^ S) ] ss 0. These are the analytical con¬ 
ditions for EM cloaking to occur in both cases (i) and 
(ii). Notice that, for this set of cloaking conditions, one 
can readily obtain case (i) from case (ii) by replacing A n 
with B n , and vice versa. In the following, we present 
a special configuration where multiple electric and mag¬ 
netic Fano resonances in the scattering cross-section are 
achieved. 
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B. Core in the Rayleigh limit (|mi|a; -C 1, x <C 1) we obtain: 


Let us consider now a non-dispersive dielectric 
core (ei,/xi) coated with a dispersive magnetic shell 
[£ 2 ( 0 ;), P 2 (u)\, with |/Z 2 r(w)| 1, both of them lossless 

media. Fixing the size parameters x = ka and y = kb, the 
cloaking condition at resonance ui = uiq [Qsca(^o) = 0] 
takes place for a certain refractive index m 2 = m 2 (^ 0 ) 


~(a,b) 


and 


tan 77722 / + A^ = cot m 2 y + A^ = 0 . (18) 


Since the arguments of these trigonometric functions are 
very large, it is convenient to consider a small frequency 
variation Sui in the vicinity of uio, leading to a corre¬ 
sponding variation Sm (which is not necessarily small) 
in the refractive index: 7712(^0 + Sui) = m 2 + Sm. To 
simplify our analysis, consider that only £2 = £ 2 (ui) 
depends on the frequency in this range Sui, so that 
the impedance index remains approximately unchanged: 
to 2 (w 0 + Scu) ft m 2 (u>o) =m 2 «l. 

In addition to the small-particle (x < y <C 1) and ferro¬ 
magnetic (\m 2 \y 3> 1 ) limits for the shell, we impose that 
the dielectric core satisfies the Rayleigh limit: \m\\x -C 
1. This provides A\ ft [£ lr a;tan(m 2 a:) + 2777 - 2 ]/[£ir£ — 
2 fh ,2 tan(m 22 ;)], where we have used Di(m\x) ft 
2/(m\x) + 0(x) and B\ is readily obtained from A\ 
by replacing (£i r , 777,2) with ( pi r , I/TO 2 ). Consequently, 

A^’^^wo + Sui) ft A^ ’ * — ( 8m)x . From the latter ap¬ 
proximation, we can finally write the analytical expres¬ 
sions in terms of the aspect ratio S = a/b = x/y, which 
is a geometric parameter that plays a crucial role on the 
EM cloaking. Imposing S ft 1, we obtain 


tan 


m 2 (ui 0 +Suj)y + A^(w 0 + Sui) 
ft tan 


(m 2 + Sm) y + A 1 — (Sm) x 


= tan 


m 2 y + A^ + 7/(1 — S)5m 


tan 

-r-(A)l 

77727/ + Al 

+ tan [7/(1 — S)Sm] 

1 — tan 

77727/ + A 

(A)] 

tan [2/(1 — S , )dm] 


= tan [//(l — S)Sm] 
ft 2/(1 — S)Sm ; 


_(.B) 

similarly, cot[m 2 ?/ + A x + 7/(1 — S)Sm ] = — tan[ 2 /(l — 
S)Sm] ft ~2/(l — S)Sm , since we can choose \Sm\y(l — 
S) < 1 even if \Sm\y 1 (small shell thickness). Hence, 
by assuming a shell with high refractive index so that 
| m 2 1 2 / ^ 1, the condition S' ft: 1 guarantees the validity of 
our analysis. Notice, however, if one imposes \Sm\y <C 1, 
this analysis is valid for every 0 < S < 1. In particular, 
for S ft 0 [homogeneous sphere (£ 2 ,^ 2 )], we retrieve the 
results of Ref. [25|. Substituting these approximations 
into Eqs. HD m and redefining the auxiliary functions, 


f[ a) 

g[ a) 

■p(B) 


2 mo 


j/ 2 (l - 5) 

m2 


+ Sm 


y 2 (l — S) 

,3 


— Sm 


r 

3 


l-S 


— 2 Sm 


(19) 

( 20 ) 
( 21 ) 



m 2 

l-S 


+ Sm . 


( 22 ) 


From Eqs. m and (HOD, the constructive and destruc¬ 
tive interferences in |ai| 2 occur for Sm ^ = m 2 /\y 2 (l — 

S)] and Sm^ = —2m.2/[y 2 (l — S)], respectively. Anal¬ 
ogously, constructive and destructive interferences in 
|6i| 2 occur for Sm ^ = — 7772/(1 — S) and Sm [= 
77 I 2 /[2(1 — 5)], respectively. In particular, we have an¬ 
alytically demonstrated that both | <3.1 1 2 and |6i| 2 have a 
Fano lineshape |ai| 2 cx (p +/3) 2 /(p 2 +1), where p = 
[(1 + a)Sm — aSm^ — 5m' A ^ > ]/\\ r )\' S /a(?> — a)], a = y 6 /9, 
/3 = 27 /[| 7 |y / a (3 — a)] and 7 = Sm^ — Sm^ A \ The 
|t»i | 2 profile is obtained from |ai | 2 by replacing Sm^ 

with Sm Ip . Assuming that only £2 depends on u> and 
\p 2 v\ ^ 1, we expect a Fano lineshape for both |ai| 2 and 
|6i| 2 as a function of £ 2 ( 0 ?). 

These asymmetric dipole resonances in total cross- 
sections are referred to as unconventional Fano reso¬ 
nances [HI, HI, which usually occur beyond the appli¬ 
cability of the Rayleigh approximation and are intrinsi¬ 
cally different from the conventional Fano resonance that 
arises from the interference between different multipole 
orders in a specific scattering direction j(| [28[ . Here un¬ 
conventional Fano resonances result from the interference 
of different EM modes excited inside the scatterer due 
to its high refractive index (hence, beyond the Rayleigh 
limit) with the same multipole moment 77 = 1. Moreover, 
note that the asymmetry parameter (5 has different signs 
in |ai| 2 and |6i| 2 profiles as a function of frequency [see 
the terms ztSm in Eq. m » This provides a config¬ 
uration in which the destructive interferences of the elec¬ 
tric and magnetic dipoles coincide without overlapping 
the scattering resonances. Therefore, we have “comblike” 
electric and magnetic dipole resonances for |rni|:r <C 1 , 
\m 2 \y ^ 1, and 2 / 1, with the additional condition 

5 ft 1 if \Sm\y 1. The analysis for the case |£ 2 r | 1, 

with |p, 2 r(w)| -C 1, is completely analogous. 


C. Comblike resonances in core-shell particles 

We have demonstrated so far that it is possible to ob¬ 
tain, in the long-wavelength regime (kb -C 1), multiple 
Fano, comblike resonances in the scattering cross-section 
of concentric core-shell spheres. The basic conditions 
consist of a core in the Rayleigh limit (\m\\ka -C 1, 
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ka -C 1 ) and a concentric shell with both large refrac¬ 
tive index (\m 2 \kb ^ 1 ) and aspect ratio near to unity 
(a ~ b). With a proper choice of parameters, one can 
obtain electric and magnetic Fano-like resonances with 
a scattering minimum at the same frequency (antireso¬ 
nance or Fano-dip). As we shall demonstrate, by making 
these approximations the Fano-dips of the electric and 
magnetic dipole resonances overlap without any addi¬ 
tional assumption, giving rise to the unconventional Fano 
resonances. However, the approximation S ~ 1 (small 
shell thickness) is now mandatory. 

To identify the positions in frequency where these 
comblike resonances occur, we set m 2 —> 0 for |/X 2 r| —► 
00 (ferromagnetic limit), leading to A\ ss tan(?B 2 a;) 
and B\ « — cot(m 2 x). From Eqs. (fill) and (flfil) . the 
cloaking condition becomes tan[? 7 i 2 y(l — 5)] = 0, i.e., 
m 2 y(l — S) = Nn, with N = ± 1 ,± 2 ,... Note that 
we impose S ss 1 because \m 2 \y ^>1. In partic¬ 
ular, this cloaking condition is readily obtained from 
Eq. ITU) , but not from Eq. (fl 6 l) . By showing that ss 
tan _ 1 [cot(m 2 a;)] = 7t/2 — cot^ 1 [cot(m 2 x)] = tt/2 — m 2 X 
and using the trigonometric identity cot(y>i + ^ 2 ) = 
(cot Lp\ cot f 2 — 1 )/(cot <£i+cot ip 2)1 one can easily demon¬ 
strate that Eq. m also provides tan[m 2 i/(l — S')] = 0 . 
Explicitly, we have: 


cot 


777 - 22 / + &[ B) 


cot [ 77122 / + tt /2 — 7772a:] 

cot [ 77722/(1 — S) + 7t/ 2] 
— tan [ 77722/(1 — S')] . 


The condition N = 0 (i.e., S = 1 since \m 2 \y 1) leads 

to a homogeneous dielectric sphere (ei, /zi), and, thereby, 
there is no cloaking. These integer numbers N, which dis¬ 
cretize EM cloaking for a certain set of parameters, can 
roughly be interpreted as indicating the “cavity modes” 
of the magnetic shell. In fact, being A 2 the wavelength in¬ 
side a nanoshell with thickness d= b—a , the EM cloaking 
occurs for A 2 « 2 d/N. Note that this is the same condi¬ 
tion for resonant modes in a Fabry-Perot cavity with d 
being the separation distance between its mirrors. Since 
this analytical result is typical of resonant optical cavi¬ 
ties, it suggests that the cloaking condition does not de¬ 
pend on the geometry of the particle (optical cavity), but 
rather the magnetic material parameters, shell thickness, 
and irradiation schemes. Analogously, similar results are 
obtained for the reciprocal configuration (high permittiv¬ 
ity coatings) setting 1 / 777-2 — > 0, for |e2r| —> 00 , where the 
same cloaking condition is achieved. 

Rewriting the cloaking condition, we explicitly have 
a mutual dependence between £ 2 r and /z 2 r , for y -C 1, 

1777 -1 \x -C 1 , \m 2 \y > 1, and S w 1: 


£2r/72r 


7 TN 

. 2/(1 — S )_ 


N = ±1, ±2,... 


(23) 


Note from Eq. (1231) that for y sa 1, large values of /i 2 r 
(or £ 2 r) are obtained for S ~ 1 (small shell thickness). 
Besides, Eq. © does not depend on (£i,/zi), and this 


independence on the core material parameters is valid 
for |mi|a; -C 1. This is of great interest, for instance, to 
cloak an arbitrary dielectric particle. Also, due to the 
electric and magnetic Fano-like resonances, the cloaking 
frequencies given by Eq. (1231) are spectrally near to a 
pair of scattering resonance peaks. Each of these scat¬ 
tering, Fano-comb resonances are associated with a mode 
number N and are referred to as morphology-dependent 
resonances of first order (n = 1). 


IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

In the following we do not restrict ourselves to the 
small-particle limit, in which our previous analytical re¬ 
sults are based, but we rather consider the exact expres¬ 
sions within the Lorenz-Mie theory, given by Eqs. m 0 . 
In particular, we consider a dielectric core with optical 
properties (£i r = 10,/zi r = 1). For the shell material, 
we consider a lossy Drude model for the relative per¬ 
mittivity, £Drude = £oo — /\u}(u + *F)], with the same 

parameters as used in Ref. [Ill] for aluminum-doped zinc 
oxide semiconductors: £00 = 3.3, w p = 2213.2 THz, and 
r = 0.002w p . The core radius is a = 100 nm and the 
aspect ratio is S = 0.9 (b « 0.13A). For the shell perme¬ 
ability, instead of assuming /Z 2 r = 1 (nonmagnetic mate¬ 
rial), we investigate the cases /z 2r = 10 3 , 2 x 10 3 , 5 x 10 3 , 
and 10 4 . 

Figure [2] shows comblike, multiple Fano resonances in 
the scattering efficiency Q sca for Re[£ 2 r (w)] > 0 and lo in 
the vicinity of w p , with 1.2 < m\x < 1.7 and 0.4 < y < 1. 
In Fig.[2][a), we plot the scattering efficiency for y .21 = 10 3 
(i.e., 1.7 < \m 2 \y < 48) showing both the electric and 
magnetic dipole contributions with the corresponding an¬ 
tiresonance or Fano-dip in between. The inset of Fig.[2](a) 
shows the scattering efficiency Qi ca P ^ calculated from the 
electric and magnetic dipole approximations, given by 
Eqs. (fTTl) 1171) . As can it be confirmed by comparing 
the inset to the main plot of Fig.(2]]a), these approxima¬ 
tions are in good agreement with the exact Lorenz-Mie 
theory, from which the main plot in Fig. [2)]a) is calcu¬ 
lated. However, from now on, all the calculations are 
performed using the exact expressions for the Lorenz- 
Mie coefficients [26]. Increasing /Z 2 r up to 2 x 10 3 (i.e., 
2.4 < 1 7772 1 y < 67), a second scattering dip appears in the 
same spectral frequency range, as evinced in Fig. [2|b). 
A very good estimate of the frequency position of these 
scattering minima is given by Eq. G3D, even though the 
Rayleigh and the small-particle limits are not entirely 
satisfied. For this set of parameters, very efficient cloak¬ 
ing is achieved for a dielectric sphere (£ ir = 10, fi\ T = 1) 
around w « 1.08u; p , for N = 2 and /Z 2 r = 2 x 10 3 . We 
emphasize, nevertheless, that our aim is not to find an op¬ 
timal set of parameters to cloak a dielectric particle [ 20 ] , 
but rather to achieve Fano-comb scattering resonances as 
a function of the magnetic permeability of the coating. 

In Figs.[3Ka) and[3]]b), we show the comblike scattering 
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Figure 2. (Color online) Scattering efficiency Q sca for a dielec¬ 
tric sphere (ei r = 10, fi i r = 1) coated with a plasmonic, lossy 
magnetic shell [e 2 r(^) = £Drude,P 2 r 3> 1]. The inner sphere 
has radius a = 100 nm and aspect ratio S = a/b = 0.9. 
The dotted line represents the sphere without coating. The 
cloaking occurs at £21 values given by Eq. (HD- (a) Q sea for 
P 2 r = 10 3 . The |cti|“ and |for| 2 curves refer to the electric 
and magnetic dipole contributions within the exact Lorenz- 
Mie theory, respectively. The inset is calculated from the 
dipole approximations (ED (HD, showing that Qsca P ' ) is in 
good agreement with Q sca . (b) Q sca for /x 2r = 2 x 10 3 . 


response for y, 2 r = 5x 10 3 and compare it to the analyt¬ 
ical prediction given by Eq. (1231) . which determines the 
values of £21 for which these multiple Fano resonances 
occur. As it can be seen from Fig. [31 Eq. (PHD gives an 
excellent prediction for the frequency position of comblike 
resonances even for a lossy shell with £ 2 r(w) > 1, y ~ 1, 
and for higher orders of N. Figure [31(b) corroborates 
that Eq. (1231) is valid even beyond the Rayleigh limit for 
the core and is robust against absorption, provided that 
\rri 2 \y ^ 1 (with I/ 2 . 2 1 ^ |^ 2 1) and S ss 1. In particu- 
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Figure 3. (Color online) Comblike scattering response for a 
dielectric sphere (ei r = 10, pi r = 1) coated with a plasmonic, 
lossy magnetic shell [£ 2 t(oj) = £Drude,/ 22 r 3> 1]- The inner 
sphere has radius a = 100 nm and aspect ratio S = a/b = 
0.9. The cloaking occurs at £ 2 r values given by Eq. P5D. 
(a) Scattering efficiency Q sca and absorption efficiency Q a b s 
(inset) for /X 2 r = 5 x 10 3 . (b) Shell permittivities £ 2 r calculated 
from the Drude model and from Eq. (1231) for p 2 r = 5 x 10 3 
and N = 1, 2, and 3. The intersection between the vertical 
lines and the Drude curve corresponds to the localization, in 
frequency, of the Fano-comb. 


lar, note that each value of N replicates the same scat¬ 
tering profile of electric and magnetic dipole resonances, 
but with different amplitudes, and that they are approx¬ 
imately equally spaced in frequency. Also, the inset in 
Fig. [3Pa) shows high absorption in these antiresonance 
scattering points (Fano-dips). This high absorption can 
be associated with resonances of the EM energy within 
the scatterer [27i,[3Q|]. I n fact; according to Refs. jil9l.f28l]. 
the Fano resonance in the scattering cross sections may 
lead to off-resonance field enhancement within the parti¬ 
cle. 

It is important to mention that the profile and the po¬ 
sition of the comblike resonances can be externally tuned 
in engineered magnetic materials. Indeed, by applying an 
external magnetic field and/or varying the temperature, 
one can change the value of y, 2 - Ferromagnetic mate¬ 
rials, for instance, typically follow the Curie-Weiss law 
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above their critical temperature T c [3lj |: /i 2 r oc (T—T c )~ l . 
Metamaterials made of ferroelectric compounds are ex¬ 
amples of systems with tunable magnetic permeability at 
terahertz frequencies, with relatively low losses [32j]. In 
these materials the tunability is achieved by varying the 
temperature, so that they could possibly be employed in 
an experimental verification of our findings. This possi¬ 
bility of tuning the value of /z 2 has influence not only on 
the resonance position but also on the number of reso¬ 
nances that compose the comblike scattering response at 
a given spectral range. This can be seen, e.g., in Fig. [2] 


We emphasize that, to the best of our knowledge, this 
tunability is not found in previous proposals involving 
comblikc resonances in core-shell nanoparticles, which 
are based mainly on nonmagnetic materials [ll,|l2|. In 
these previous systems, the position and profile of comb¬ 
like resonances are determined a priori by the nanostruc¬ 
ture design, either by varying the plasma frequencies with 
a gradient of doping level in multilayered spheres [ll| or 
by constructing nonconcentric core-shell nanowires [121 ]. 
This fact, which distinguishes our proposal from the pre¬ 
vious ones, allows for tunable, versatile comblike scat¬ 
tering response in concentric, single-layered core-shell 
nanoparticles. 


In Fig. Ha), the shell permeability /z, 2 r = 10 4 is one 
order of magnitude greater than in Fig. Ha), which im¬ 
plies that more scattering dips occur at the same spectral 
range. The differential efficiencies for the back and for¬ 
ward scattering [ 2 ^], respectively, Qback = | S^Li( 2n + 
l)(-l)”(a n - b n )\ 2 /y 2 and Qforward = IE^i( 2rl + 
1 )(a n + b n )\ 2 /y 2 , are plotted in Fig. Ht>). Due to the 
fact that the shell is magnetic (/J 2 r 7 ^ 1 ) light scattering 
is asymmetric with a backscattering dominance; indeed, 
Qback has essentially the same profile as < 3 sca , as can be 
verified in Figs. H a ) and HT*)■ On one hand, the inter¬ 
ference between electric and magnetic dipoles facilitates 
the suppression at the scattering dips of order N in the 
backscattering direction. On the other hand, the for¬ 
ward scattering efficiency Qforward, plotted in the inset 
of Fig. Hb); shows a strong suppression in the spectral 
frequency region between the N and TV + 1 dips, and 
additional Fano-dips due to the conventional Fano res¬ 
onance. This strong asymmetry between back and for¬ 
ward scattering may be exploited in applications. In¬ 
deed, the ultrasharp Fano-comb scattering suggests ap¬ 
plications in optical tagging and signal processing, for in¬ 
stance. In contrast to previous studies on multiple Fano 
resonances [ll, ]T 2 j] , our approach not only allows one to 
analytically predict, by Eq. (1231) . the frequency position 
where these resonances are expected to occur, but also 
allows for the possibility of tuning their positions via the 
external variation of the magnetic permeability of the 
nanoshell. 
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Figure 4. (Color online) Comblike scattering response for a 
dielectric sphere (ei r = 10, /in = 1) coated with a plasmonic, 
lossy magnetic shell [e 2 r(w) = £Drude,/t 2 r = 10 4 ]. The inner 
sphere has radius a = 100 nm and aspect ratio S = a/b = 0.9. 
The cloaking occurs at £21 values given by Eq. (1231) . (a) Scat¬ 
tering efficiency Q sca . (b) Differential scattering efficiencies 
Qback and Toward (inset). 


V. CONCLUSIONS 

Within the Lorenz-Mie theory, we have investigated 
the possibility of achieving multiple Fano resonances in 
the scattering response of dielectric spheres coated with 
a single-layer, concentric magnetic shell. We have de¬ 
rived, in the long-wavelength limit, explicit analytical 
conditions for the occurrence of comblike resonances in 
core-shell spheres. In particular, we have shown that 
nanoshells with high magnetic permeability or high elec¬ 
tric permittivity values induce the formation of multiple 
Fano, comblike resonances composed of electric and mag¬ 
netic dipole resonances. These Fano-comb resonances 





















have been shown to follow a typical relation of a reso¬ 
nant optical cavity. As the condition for the occurrence of 
multiple Fano, morphology-dependent resonances explic¬ 
itly depends on the magnetic permeability of the shell, 
we argue that both the position and profile of the comb¬ 
like resonances could be tailored by an external mag¬ 
netic field and/or by varying the temperature. Together 
with the simplicity of employing a single-layered core¬ 
shell nanoparticle, this tunability of the single scatter¬ 
ing response distinguish our proposal from the previous 
ones considered so far, which are based on e.g. multilay¬ 
ered semiconductor materials or nonconcentric core-shell 
particles. These singular scattering properties make the 
system proposed here a tunable, versatile optical device 
that may find applications in multifrequency biosensing, 
optical tagging and signal processing. 
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